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1 Introduction 



The cascade model (Cohen (1990), Cohen and Newman (1985, 1986), Cohen, 
Briand and Newman (1990), Newman (1992)) is introduced to explain the eco- 
logical data on community food webs. The model generates a food web at 
random. In it the species are labeled from to m, and arcs are given at random 
between pairs of the species. For an arc with endpoints i and j (i < j), the 
species i is eaten by the species labeled j. The chain length (height) is compared 
with the length of food chain in ecological data. The problem of chain length 
of food webs gives natural questions for random evolution of graphs in which 
connectedness and other related problems are studied (Erdds and Renyi (1960), 
Shepp (1989), Durrett and Kesten (1990)). 

The random sequential bisection of intervals (Sibuya and Itoh (1987), Jan- 
son and Neininger (2008), Dutour Sikiric and Itoh (2011) ) gives an analogous 
asymptotic behavior to the binary search tree (Robson (1979), Flajolet and 
Odlyzko (1982), Mahmoud and Pittel (1984), Devroye (1986)). Here we in- 
troduce a random sequential generation of intervals to understand the chain 
length for the cascade model, extending the idea of random sequential bisection 
of intervals. We can naturally obtain a formula for the chain length, as given 
in section 2 and section 3, which suggests the corresponding formula for the 
original cascade model as given in section 4. The random sequential interval 
generation is a Poisson approximation for the cascade model of food webs, in 
the sense given in equation (|12p in section 3, and equation (|16l) in section 4. 

Consider the random oriented graph with vertex set {0, 1, 2, m} in which 
the ( m ^ _1 ) oriented edges (j, j) with i < j occur independently of each other with 
probability p — c and no edge j < i occurs. Let us give a random variable X^ 
to each oriented edge with i < j with the vertex set {0, 1, 2, m}, where 
arc mutually independent random variable with X^ = 1 with probability 
p = c and X^ = with probability 1 — c. 

Consider the random variables Xij with i = 0. For sufficiently small c, the 
number N = Xij with i = is approximately distributed by the Poisson 

distribution ^-^—e~ cm . Let each of N random variables -X"o,fc(i)> ^o,fc(2)i — , Xo t k{N) 
take the value 1. Consider for each k(l), I = 1,2, ...,N , the random variables 
Xk(l),k(l)+i, -Xfe(j),fc(J)+2> ...,-Xfc(i),m- Then Y^jLk(i)+i x ij witn * = ^(0 is a P~ 
proximately distributed by the Poisson distribution ( c ( m ~"^( ; )) e -c(m-fc(0) < w e 
define the Poisson generation of intervals to define the random sequential gen- 
eration of intervals for the cascade model as follows. 

Poisson generation of intervals is defined for the interval [0, y] to take a 
random variable N(y) distributed by the Poisson distribution with the param- 
eter cy as 

P(N(y) = k) = ±(c y ) k e- c y 

to generate the intervals [0,Xi(y)], [0, X2(y)], [0, XN( y )(y)], where Xi(y) is 
distributed uniformly at random on the interval [0,y], mutually independently, 
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for each i. 

At step 1 we apply the Poisson generation to the interval [0, a;]. At step 
1 < j, for each interval [0,y] generated at the step j — 1 apply the Poisson 
generation to the interval independently from other intervals and independently 
from the previously generated intervals. Each interval which does not generate 
any interval at step j does not generate any interval after the step j. We continue 
the steps as long as we have at least one generated interval. 

We also define the random sequential generation of intervals with an expo- 
nentially distributed starting point as follows. At step we generate the interval 
[0, x — Z] where Z is distributed by the density e~ cz . If ([0, x — Z] <£_ [0, x}), we 
finish and stop to make the next step. If ([0, x — Z] C [0, x]), we proceed to step 
1. At step 1 we apply the Poisson generation to the interval [0, x — Z]. At step 
1 < j, for each interval generated at the step j — 1 apply the Poisson generation 
to the interval independently from other intervals and independently from the 
previously generated intervals. Each interval which does not generate any in- 
terval at step j does not generate any interval after the step j. We continue the 
steps as long as we have at least one generated interval. In the original cascade 
model of food webs, the value cm is assumed to be a constant. Assuming c is 
a constant independent from m, the asymptotic behavior of the longest chain 
is of mathematical interest, which may have applications for example to task 
graphs for parallel processing in computer science (Newman (1992)). As shown 
in sections 2, 3, 4, and 5, the random sequential generation of intervals helps to 
understand the asymptotic length of the longest chain of the cascade model. 



2 Random sequential generation by the Poisson 
distribution 

For each stopped interval generated by the above procedure, we count the num- 
ber of steps to get the stopped interval. Let K(x, a) be the number of intervals 
which take a steps until stopping. Let L(x,a) be the expectation of K(x,a). 
We get the asymptotic behavior of L(x, a) as in the case of random sequential 
bisection for the binary search tree (Sibuya and Itoh (1987)). 
We have 

L(x,0) = e~ cx . (1) 

For 1 < a, we have, 

I r x 00 I 

L(x,a) = - ^-{cx) k e- cx kL(y,a-l)dy, 

X Jo k=1 K. 

= - / ET i (o)'e-°fct(l/,a-l)ii!l ) 
x Jo k=1 «• 
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= a;- / V-( ra ) fc e - cx L(y, ffl - l)dj/, 

- c / Lfo.a-l) dy. (2) 
Jo 



Hence for a — 1, we have 



= c / L(y,0) dy, 
Jo 



= c [ e~ cy dy. 
Jo 



We have finally 



px rXa-l rxi 

L(x,a) — c c ...c / c~ CXl dx\ dx 2 ...dx a -i 

Jo Jo Jo 



Hence we have 



(rr) a 

L(x,a) + L(x,a + 1) = ^-. (4) 



Put a — kx, and apply the Stirling formula n! ~ v / 27m™ +1//2 e ™, then 

L(x,a) + L(x,a + l) = ^j^~ _ ^f^* — . (5) 

(fcx)! V^r(fcz) fex+1 / 2 e- fe * v 7 



('■■'•)*■' ! ; (£)^ e ^ ( 6 ) 



As 



2ir(kx) kx+1 / 2 e- kx y/2irkx k' 

we have the following theorem, which corresponds to the stronger result by 
Newman (1992) for the cascade model. 

Theorem 1 

(i) 

(rr) a 

L(x,a) + L(x,a + 1) = ^-. (7) 

(ii) Put a = k x, for x — > 00, if e c < k, 

L(x,a) +L(x,a + 1) ->• 0, (8) 

if fc < e c, 

L(x,a) + L(x,a + 1) -> 00. (9) 
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3 The chain length with an exponentially dis- 
tributed starting point 

The expected chain length for the case that the starting interval is [0, x — u] 
(c [0, a;]) where u is distributed by the density ce~ cu is given by 

r 

M{x,a) = I L{x — u,a)ce cu du. 
Jo 



M(x,a) + M(x,a+l) = \ (L(x - u, a) + L(x - u, a + l))ce~ cu du 

v~ du , (10) 



M(x,a)+M(x,a + 1) = (-1)^* + ]T(-1)^ ^4 (11) 

Hence we have Theorem 2 in the same way to Theorem 1 by using the Stirling 
formula. 

Theorem 2 

(i) 

M(x,a)+ 2M(x,a + l)+M(x,a + 2) = ^ ^ . (12) 

(ii) Put a = k x, for x oo, if e c < k, 

M(x, a) + 2 M (x, a + 1) + M (x, a + 2) 0, (13) 

if k < ec, 

M (x, a) + 2 M(x, a + 1) + M (x, a + 2) -> oo. (14) 



4 On Cascade model 

Putting m = 5— 1, the expected number of chains with length n on the cascade 
model (Cohen and Newman (1986)) is given by 

E(C n )= P -J2(S-k)( k -% s - k -\ (15) 



k—n 



By using this formula, we can extend equation (1121) to the cascade model. The 
following theorem will show that our Poisson interval tree is a natural continuous 
model for the cascade model of food webs and gives an elementary approach to 
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the theorems for the cascade model (Newman (1992)) on the expected length 
of a chain and the asymptotic length of the longest chain. 



Theorem 3. 

(i) 



E( C'n I - 2 E[ ( ',,-. J ) + EiC„-- 2 ) : : \ ^ ^ l>" C lCi ) 



(ii) Put n — k S p, for S p — > oo, if e < A;, 

£7(C„) + 2£(C n+ i) + £(C*„ +2 ) -> 0, (17) 

if k < e , 

£(C„) + 2 £(C n +i) + £(C„+ 2 ) -)■ oo. (18) 

Proof. By using 

/ C n . v ^ , ~ , , / k 1 



k—n 



^)=E(5-*)(;:ij9 fl -*- 1 . as) 



S-l 



= e a-^-^r 1 

^ k=n+l v 



gfiT-fc-1 



^-l 



E ( 5 - fc ) 



fc=n+l 
S-l 



??■ 



£ V-*, (20) 



and 



^(%) - E (i-qns-k)( k -] 

p n z — ' \ n + 1 



fc=n+2 
S-l 



q S-k-l 



E . I 

\n + 1 



fc=n+2 
S-l 



jfe-1 



k - 1 



^s fc 



fc=n+2 v 

+ e ^- fc )(";i)^ fe+1 ' (2i) 

fc=n+2 V / 

we can easily obtain eq. (fl6|) . 

(ii) is obtained by using the Stirling formula as in the proof of Theorem 1 
and Theorem 2. 
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5 On the height of the generated tree 



For each stopped interval generated by the above procedure in the section, we 
count the number of steps a to get the stopped interval. Consider the maximum 
H (x) of the numbers and let us call it the height of the generated tree. For the 
probability F(x, ft) = Pr(H(x) < ft) of the height of Poisson interval tree H(x), 
we have for h = 0, 

F(x,h)=e- cx , (22) 

for 1 < ft, 



F(x,h) 



= e~ cx 



+ Y.-TT e ~ CX — •••/ F(yi,h-l)---F(y k ,h-l)d yi ---dy k , 
k=1 K - x Jo Jo 

^-(c/ F{y,h-l)dy) k 

t-n K - Jo 



Jc f* F(y,h-l)dy) 



(23) 



We can sequentially integrate at each step and obtain F(x, h) starting from 
F(x,0) = er cx . The expansion at x = is obtained by Mathematica up to 
(cx) for example for h = 0,1, 2, 7, which gives reasonable numerical values 
otF(x, h)—F(x, ft— 1) when cx < 2.1, while it does not give reasonable numerical 
values for 2.2 < cx, For 4 < ft the first five terms of the probability F(x, ft) — 
F(x,h— 1) that the the height is ft is given by 

Jcx) h+1 (Cx) h +* (CX)^ (cxf+± 

ft! (ft + 1)! (ft + 2)! (ft + 3)! (ft + 4)!' [ ' 
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